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| Abstract. The effect of microwave irradiation on the phase dynamics of intrinsic Josephson 

r—{ _ junctions in high temperature superconductors is investigated. We compare the current-voltage 

r/3 ■ characteristics for a stack of coupled Josephson junctions under external irradiation calculated 

in the framework of CCJJ and CCJJ+DC models. 



1. Introduction 

Since the discovery of the intrinsic Josephson effect in high temperature superconductors (HTSC) 
such as BiiSr2CaCu20 y (BSCCO), numerous studies of the basis nature, nonlinear dynamics 
and device applications of the superconductors have been extensively developed. The high- 
frequency response of Intrinsic Josephson junctions (IJJs) have been constantly updated. One of 
the most spectacular indications of the Josephson effect is locking of Josephson oscillations to the 
frequency of external microwave (MW) irradiation. As the Josephson frequency is proportional 
to voltage, the locking leads to appearance of steps at quantized voltages, called Shapiro steps 
[HE]. The steps appear in the current- voltage characteristic (CVC) at discrete voltage values 
V n = nhu/2e (u, the frequency of the applied MW signal; h, the Planck constant; e, the 
elementary charge; n, integer number). 

It was shown by Koyama and Tachiki [3j that the system of equations for capacitively coupled 
Josephson junctions has a solution corresponding to the longitudinal plasma wave (LPW) 
propagating along the c-axis. So, the Josephson oscillations can excite the LPW by their 
periodical actions. The frequency of Josephson oscillations uij is determined by the voltage 
value in the junction, and at loj = 2ulpWi where lolpw is LPW frequency, the parametric 
resonance (PR) is realized. External radiation essentially changes the physical picture of the 
coupled JJs. As known, the one-dimensional models with coupling between junctions capture 
the main features of real IJJs, like hysteresis and branching of the CVC, and help to understand 
their physics. So it is of interest to investigate this system under MW irradiation using these 
models. 

In this article we investigate the phase dynamics of the coupled Josephson junctions under 
MW irradiation in the framework of both CCJJ and CCJJ+DC models. 



2. Model and Method of Calculations 

Because of the thickness of the S-layers in HTSC is extremely small and comparable to the 
Debye length rjj of the electric charge screening, therefore, there is no complete screening in a 
separate S-layer. It leads to the generalized scalar potential <£/ of S-layer. The &i is expressed 
through the electric scalar potential <f>\ and the derivative of order parameter's phase Q\ of S-layer 
I by $,(t) =4>i- Vfi. 

Without screening of the charge in S-layer the Josephson equation is generalized ipi = 
Vi + a(2Vi — Vi-i — Vi+i) (where a is the coupling between junctions) and together with the 
expression for total current through the stack I = Vi + sin ipi + (3V[ describe the phase dynamics 
of the coupled system of JJs in the framework of CCJJ model [3j. In CCJJ-model, the system 
of equations for the gauge-invariant phase differences <pi{r) between S-layers I and I + 1 in the 
presence of MW irradiation is described by: 



dP dip 

-j^tpi = (I ~ sin ipi - P~rp) + a(sin^+i + sin<^_i - 2sinc^) - As'mujT (1) 

where / is the external bias current, A is radiation amplitude, oj is MW radiation frequency and 
j3 is the dissipation parameter (/3 related to the McCumber parameter j3 c by /3 C = 1//3 2 ). The 
CVC in CCJJ model are characterized by branching at / = I c (where I c is the Josephson critical 
current) and strong branching in the hysteresis region: the number of branches in general case 
is equal to 2 N , where N is the number of junctions in the stack. 

In CCJJ+DC model [I] the diffusion current J l D = <J?m_i — $i is added and the expression 
for current through the stack of junctions in this model is I = Vi + sin <pi + fiipi. Second 
order differential equation for <pi has additional terms in compare with CCJJ model which are 
proportional to the product of coupling and dissipation parameters 

d2 n a d{ f l \ , ( ■ , ■ o • \ 

■^2<Pl = {I ~ smipi - P-^r) + a(sm^ + i + sm^_i - 2sm^) 

+ap{ — 1 2— — )— ismuT (2) 

ut dr dr 

To find the CVC of the stack we solve the system of dynamical equations for phase differences 
using the fourth order Runge-Kutta method. We use a dimensionless time r = toj p , where uj p is 
the plasma frequency u p = \/2eI c /hC and C is the capacitance. In our simulations we measure 
the voltage in units of Vq = hu p / (2e), the frequency in units of u p , the current and the amplitude 
of radiation in units of I c . To study the time dependence of the charge in the S-layers, we use 
the Maxwell equation eeoV .E =Q, where e and Eq are relative dielectric and electric constants. 
The charge density Qi in the S-layer I is proportional to the difference between the voltages Vi 
and Vi + \ in the neighbor insulating layers Qi = Qoa(Vi+i — Vi), where Qq = eeoVo/rjj. 

We solve the system of dynamical equations for phase differences at fixed value of bias current 
I in some time interval (0, T m ) with the time step 5t. This interval is used for time averaging 
procedure. Then we change the bias current by 51, and repeat the same procedure for the 
current I + 51 in new time interval (T m , 2T m ). The values of the phase and its time derivative 
in the end of the first time interval are used as initial conditions for second time interval and 
so on. The total recorded time is calculated as t r = tuj p + T m (Io — I)/ 51, where Iq is an initial 
value of the bias current for time dependence recording. In our simulations we put T m = 1000, 
5t = 0.05, 51 = 0.0001 and we add to the bias current a small noise with amplitude in the 
interval (+10 -8 , — 10 -8 ). The details concerning the numerical procedure are given in Ref. [S]. 



3. Results and Discussion 

In our present work, we describe the case w > 2ulpw-, when the SS is above PRR in CVC. It 
was shown in Ref. [6] that an increase of the amplitude of irradiation leads to the appearance 



of an additional PR above SS called as radiation related parametric resonance (rrPR). 

In Fig. (U^a) we show the CVC of the stack with 10 coupled JJs under microwave irradiation 
with u> = 2 and A = 0.5 simulated in the framework of both the CCJJ and CCJJ+DC models 
at a = 0.05, /3 = 0.2 and periodic boundary conditions. At w = 0, PR is characterized 
by breakpoint current Ib p — 0.28, and breakpoint voltage Vb p ~ 11.5, corresponding to the 
Josephson frequency ojj = 1.15. 

As we can see, the first SS is developed on the outermost branch of CVC at V = ujj * N = 20. 
Dashed line stresses this fact. Both models are coincide with the increase in bias current. 
The difference between CCJJ and CCJJ+DC models appears in the branching part of CVC in 
hysteretic region. The CCJJ model demonstrates more intensive branching and larger hysteresis 
region in compare with CCJJ+DC model. This feature could be explained by the role of diffusion 
current [7]. 




Figure 1. CVC of the stack with 10 coupled JJs under microwave irradiation with u) = 2 and 
A = 0.5 simulated in the framework of the CCJJ and CCJJ+DC models at a = 0.05, (3 = 0.2 
and periodic boundary conditions; (b) The enlarged part of the rrPR above SS on the outermost 
branch of CVC. 

As we mentioned previously, an increase of the amplitude of irradiation leads to the 
appearance of rrPR. Fig. QJb) shows (for both models) an enlarged view of the rrPR region 
on the outermost branch of the CVC (the temporal oscillations of the charge in this region are 
shown in Fig. [2]). We see that rrPR is started at higher currents for CCJJ than CCJJ+DC 
model. In addition, we observe an interesting structure before SS in CCJJ-model. 

Fig. [2ja) demonstrates the temporal oscillations of the charge in the growing region of rrPR 
for the first S-layer of the stack with ten coupled JJs at a = 0.05, (3 = 0.2, u> = 2 and A = 0.5, 
combined with CVC of the outermost branch, calculated in the framework of CCJJ+DC model. 
The filled squares mark the bias current steps in CVC. We see in this figure the correlation 
between the charge-time dependence and CVC. The irradiation change the character of the 
charge-time dependence essentially and bring about a "bump" structure on the outermost branch 
of CVC, as shown in the figure. Analogies features in CVC were recently observed experimentally 
[8]. The inset shows the charge distribution among the layers in the growing region. The charge 
on the neighbor layers is equal in magnitude and opposite in sign. This is true for all adjacent 
layers and corresponds to the 7r-mode i. e. the wavelength of the created LPW is equal to 
A = 2d. 

The time dependence of the charge in the rrPR region calculated in the framework of CCJJ- 
model is presented in Fig. [2(b). We see that the charge on the S-layer appears at a higher 



value of current (Ib p = 0.4809) in comparison with the current value in case of CCJJ+DC-model 
(hp = 0.4739). The inset shows the charge distribution among the layers in the growing region, 
illustrates that the wavelength of LPW is the same as for CCJJ+DC-model (A = 2d). 




Figure 2. Charge time-dependence in the radiation related parametric resonance region of the 
stack with ten coupled JJs at a = 0.05, f3 = 0.2, u = 2 and A = 0.5 calculated in the framework 
of (a) CCJJ+DC-model; (b) CCJJ-model. The insets show the charge distribution among the 
S-layers in the growing region. 



4. Summary 

We presented the CVC of the stack with 10 coupled Josephson junctions under microwave 
irradiation calculated in the framework of CCJJ and CCJJ+DC models. The CVC demonstrate 
different features in the radiation related parametric resonance region above SS. These results 
stress the essential role of the diffusion current in formation of the CVC of intrinsic Josephson 
junctions in high temperature superconductors. We showed that the radiation might bring to 
the "bump" structures in CVCs recently observed experimentally. 
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